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摘 要 : 针对 小 麦 植株 分 坦 多 、 器 官 间 交叉 遮挡 严重 ， 难 以 用 图 像 或 点 云 准确 提取 植株 和 器 官 表 型 的 问题 ， 
本 研究 提出 了 基于 三 维 数字 化 的 小 麦 植株 表 型 参数 提取 方法 。 首 先 提出 了 小 麦 植株 各 器 官 数字 化 表达 方法 ， 
制定 了 适用 于 小 麦 全 生育 期 的 三 维 数字 化 数据 获取 规范 ， 并 依据 该 规范 进行 数据 获取 。 根 据 三 维 数字 化 数 
据 的 空间 位 置 语义 信息 和 表 型 参数 的 定义 ， 提 出 了 小 麦 植株 表 型 参数 计算 方法 ， 实 现 了 小 麦 植株 和 器 官 长 
度 、 粗 度 和 角度 等 3 类 共 11 个 常规 可 测 表 型 参数 的 计算 。 进 一 步 提 出 了 定量 描述 小 麦 株 型 和 叶 形 的 表 型 指 
标 。 其 中 ,植株 围 度 通过 基于 最 小 二 乘法 拟 合 三 维 离散 坐标 计算 ,用 于 定量 化 描述 小 麦 植株 松散 /紧凑 程度 ; 
小 麦 叶片 卷曲 和 扭曲 程度 为 定量 化 叶 形 的 指标 ， 根 据 叶 面 向 量 方向 变化 计算 得 到 。 利 用 丰 抗 13 号、 西 农 979 
号 和 济 麦 44 号 三 个 品种 小 麦 起 身 期 、 拔 节 期 、 抽 穗 期 三 个 时 期 的 人 工 测量 值 和 提取 值 进行 验证 。 结 果 表 明 ， 
在 保持 植株 原始 三 维 形态 结构 的 前 提 下 ， 提 取 的 茎 长 、 叶 长 、 茎 粗 、 茎 叶 夹 角 与 实测 数据 精度 相对 较 高 ，R 
分 别 为 0.93、0.98、0.93、0.85; 叶 宽 和 叶 倾 角 与 实测 数据 的 记分 别 为 0.75、0.73。 本 方法 能 便捷 、 精 确 地 提 
取 小 麦 植株 和 器 官 形 态 结构 表 型 参数 ， 为 小 麦 表 型 相关 研究 提供 了 有 效 技术 支撑 。 
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一 、 


1 518 是 植物 功能 -结构 模型 “以 及 数字 植物 技术 等 


A 
m 研究 的 核心 内 容 。 随 着 各 学 科 的 发 展 和 数据 处 理 


官 、 植 株 和 群体 结构 及 功能 特性 的 物理 、 生 理 和 代 "1 ， 传 统 测量 方式 已 无 法 满足 植物 表 型 高 通 
生化 性 状 ， 其 本 质 是 植物 基因 图 谱 的 时 序 三 维 表 量 获 取 和 需求 中， 迫切 需要 针对 不 同 植物 开展 精 
达 及 其 地 域 分 异 特征 和 代 际 演进 规律 "。 植 物 表 准 、 高 效 表 型 高 通 量 获取 研究。 小 麦 是 世界 三 大 


型 研究 是 遗传 育种 ”和 精准 栽培 ”的 关键 ， 也 ”粮食 作物 之 一 ， 开 展 小 麦 植 株 形态 结构 表 型 参数 
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提取 方法 研究 对 于 小 麦 株 型 育种 ”具有 重要 意义 。 

传统 的 植物 形态 结构 表 型 获取 以 人 工 测 量 为 
主 ， 所 获取 表 型 参数 精度 受 环境 和 人 为 因素 影响 
较 大 。 近 年 来 ， 基 于 图 像 的 植物 表 型 提取 方法 获 
得 广泛 应 用 “”， 但 利用 单 幅 图 像 难以 解决 植物 
带 官 复杂 的 遮挡 问题 ， 部 分 形态 结构 参数 的 提取 
精度 受 图 像 获 取 和 角度 影响 较 大 。 相 比 之 下 ， 基 于 
三 维 数据 提取 植物 形态 结构 表 型 可 以 得 到 较 高 的 
精度 “”。 在 小 麦 三 维 表 型 提取 方面 ， 付 唱 波 
等 "利用 多 视角 图 像 采集 系统 获取 小 麦 幼苗 图 
片 并 提取 了 小 麦 的 株 高 和 叶 长 ， 但 仅 针对 小 麦 分 
幕 前 的 幼苗 期 进行 了 试验 ; 翟 苗 苗 ”利用 Li- 
DAR 获取 小 麦 植株 三 维 数据 并 解析 出 高 度 和 体 
只 性 状 ， 用 于 构建 小 麦 生物 量 模 型 。 可 以 发 现 ， 
对 于 小 麦 这 种 分 喜多 、 叶 片 等 带 官 遮挡 严重 的 作 
物 ， 利 用 常规 的 二 维 图 像 和 三 维 点 云 解析 方法 都 
难以 提取 到 高 精度 的 植株 表 型 参数 ， 且 无 法 利用 
一 组 数据 提取 到 所 需 的 多 个 表 型 参数 ， 尤 其 是 小 
麦 生育 后 期 植株 内 部 的 形态 结构 信息 。 三 维 数字 
化 仪 是 一 种 可 以 测量 物体 特征 点 精确 空间 坐标 位 
置 的 电磁 设备 ， 特 别 适 用 于 分 校 结构 复杂 的 植 
物 “”。 为 解决 器 官 间 存在 大 量 交叉 逛 挡 植物 形 
态 结构 表 型 高 精度 提取 问题 ， 研 究 者 利用 三 维 数 
字 化 仪 (3D Digitizer) 直接 还 原 植物 植株 在 三 维 
空间 中 形态 结构 ， 并 在 此 基础 上 提取 植株 表 型 参 
数 。 例 如 ， 温 维 亮 等 ” 针对 玉米 形态 结构 特 
征 ， 制 定 了 玉米 三 维 数字 化 数据 获取 标准 ， 并 在 
此 基础 上 实现 了 玉米 株 型 参数 的 准确 解析 ; 
Fournier E ”利用 三 维 数字 化 仪 对 小 麦 叶片 出 生 
方向 、 叶 脉 、 叶 形 进 行 准 确 测 量 并 构建 了 数据 
集 ， 并 用 于 构建 小 麦 植株 生长 发 育 模型 。 然 而 ， 
目前 未 见 利 用 三 维 数字 化 仪 对 小 麦 植株 进行 精确 
的 形态 结构 表 型 解析 的 报道 。 

本 研究 以 此 为 出 发 点 ， 和 针对 不 同 生育 期 小 麦 
植株 形态 结构 特征 ， 利 用 三 维 数字 化 技术 ， 制 定 
小 麦 植 株 三 维 数字 化 数据 获取 规范 ， 从 “ 单 株 - 
单 茎 -器 官 ” 不 同 尺度 一 次 性 快速 提取 所 需 的 多 
个 小 麦 植株 表 型 参数 并 进行 可 视 化 。 为 更 加 深入 


分 析 小 麦 三 维 形态 结构 ， 本 研究 在 玉米 株 型 参数 
提取 方法 ” 的 基础 上 ， 针 对 小 麦 的 复杂 分 碍 结 
构 进 一 步 详 细 了 数据 获取 规范 。 利 用 该 规范 获取 
的 小 麦 三 维 数字 化 数据 相 较 于 骨架 信息 ， 包 含 了 
叶片 精细 形态 信息 ， 能 提取 更 多 三 维 表 型 参数 。 
同时 ， 还 对 植株 松散 程度 和 叶片 卷曲 、 扭 曲 程 度 
进行 了 定量 化 计算 。 在 此 基础 上 ， 进 行 小 麦 植株 
的 下 一 步 三 维 重 建 和 表 型 解析 能 够 解决 絮 官 间 和 让 
挡 严 重 的 问题 ， 从 而 获得 更 准确 、 完 整 的 信息 。 


2 小 麦 植 株 三 维 数字 化 数据 获取 与 
表 型 参数 提取 方法 


在 不 破坏 小 麦 原 始 三 维 形态 结构 的 前 提 下 ， 
对 大 田 种 植 的 小 麦 植株 进行 取样 ， 并 依据 制定 好 
的 数据 获取 规范 进行 小 麦 形态 数据 获取 。 获 取 完 
毕 后 ， 在 Python3.8 环境 下 利用 Open3D 0.13.0 对 
不 同 生育 时 期 的 小 麦 植株 进行 可 视 化 。 同 时 进行 
相关 表 型 参数 的 自动 化 提取 ， 包 括 高 度 和 长 度 参 
数 、 粗 度 参 数 、 角 度 参数 ,以 及 叶片 卷曲 和 扭曲 
程度 四 类 。 图 1 给 出 了 小 麦 表 型 参数 提取 方法 的 


整个 流程 。 


| 


数据 获取 设备 


长 度 和 高 度 参数 


VIPS 
y Bi : N 
粗 度 参数 é i 
EBB Bs 
叶片 卷曲 、 git | i 
程度 参数 “ 表 型 参数 自动 化 提取 ”数据 可 视 化 
Bl 小 麦 表 型 参数 提取 整体 流程 


Fig. 1 Overall process of wheat phenotypic 


parameter extraction 
2.1 小 麦 三 维 数字 化 数据 获取 


2.1.1 数据 获取 
冬小麦 试验 地 点 为 北京 市 农林 科学 院 试验 田 
(N39°56', E116°16)。 选 取 了 丰 抗 13 号 (FK13)、 


152 


ChinaXiv 合 作 期 刊 


202302.00157v1 


chinaXiv 


西 农 979 CXN979) 和 济 麦 44 号 (JM44) = 
个 株 型 、 叶 形 差异 较 大 的 品种 。 其 中 ， 丰 抗 13 
号 植株 较 高 ， 株 型 较 松 散 ， 叶片 具有 细 长 的 形态 
特征 ; 济 麦 44 号 植株 较 矮 ， 株 型 较 紧 凑 ， 叶 片 
具有 宽 和 短 的 特点 ; 西 农 979 号 株 型 和 叶 形 特征 
均 居中 。 播 种 时 间 为 2020 年 10 月 4 日 ， 每 个 品 
种 一 个 小 区 ， 各 小 区 大 小 设 定 为 2.25 mX 1.5 m, 
小 区 植株 行距 和 株距 分 别 为 0.2 和 0.05 mo HM] 
播种 前 撒 施 复合 肥 (N-P,O,-K,0:12-18-15)， 施 
肥 量 为 75 gm"， 并 于 拔节 期 随 水 追 施 30 g/m? PR 
Ko WATER KR, BOT. WRH., RA 
Waak, WEKEN 45 m/m. 

采用 机 械 臂 式 数 字 化 仪 Microscribe i 
(Revware，USA， 图 2) 进行 数据 获取 ， 其 有 效 
获取 范围 达 1.27 m， 精 度 为 0.178 mm， 满 足 小 麦 
植株 三 维 数字 化 数据 获取 要 求 。 


图 2 Microscribe ji 三 维 数字 化 仪 


Fig. 2 3D Digitizer Microscribe i 


于 小 麦 生长 的 三 个 关键 生育 期 (起身 期 、 拔 
节 期 和 抽穗 期 ) 对 各 品种 植株 进行 三 维 数字 化 数 
据 获取 ， 每 个 品种 重复 3 次 。 通 过 整 株 移植 至 花 
盆 、 搬 至 室内 的 方式 取样 ， 尽 量 保证 所 取样 品 为 
叶片 表面 完整 的 植株 ， 在 无 风 的 环境 下 根据 制定 
的 数据 获取 规范 对 小 麦 地 上 部 进行 三 维 数字 化 数 
据 采 集 。 数 据 获 取 前 ， 基 于 电磁 定位 原理 使 用 定 
标 块 进行 位 置 标定 ， 植 株 摆 放 在 机 械 臂 能 触 达 的 
范围 内 。 为 保持 同一 植株 各 器 官 在 同一 坐标 系 
下 ,数据 获取 过 程 中 植株 不 能 移动 。 同 时 ， 对 各 
带 官 形态 结构 参数 进行 人 工 测量 并 记录 。 在 数据 
采集 过 程 中 ， 要 求 无 风 且 仪器 不 能 触 碰 小 麦 杆 
株 ， 以 避免 小 麦 植株 叶片 在 数据 获取 过 程 中 发 生 


位 置 变 化 形成 误差 。 每 株 小 麦 植株 数据 获取 完成 
后 ， 需 利用 自主 开发 的 程序 对 数据 以 可 视 化 的 方 
式 进行 检查 ， 若 发 现 数据 缺失 或 有 误 ， 及 时 重新 
获取 。 
2.1.2 小麦 植株 数字 化 表达 及 数据 获取 规范 
基于 小 麦 生长 发 育 规 律 制定 适用 于 整个 生育 
期 的 数据 获取 规范 ， 以 确保 整个 植株 三 维 数字 化 
数据 获取 的 完整 性 和 各 结构 单元 之 间 连 接 的 正确 
性 。 草 卫星 和 李 存 东 2 利用 英文 名 称 首 字母 对 
小 麦 各 器 官 命 名 ， 系 统 地 提出 了 地 上 部 器 官 序列 
化 命名 方案 ， 且 通过 数字 标记 表明 植株 上 各 器 官 
发 育 的 时 空 序 位 。 本 研究 结合 冬小麦 的 结构 单元 
划分 ， 在 此 基础 上 以 器 官 为 单位 ， 制 定 了 小 麦 植 
株 各 器 官 数 字 化 表达 及 数据 获取 规范 。 

(1) 4. HS (Stem) 表示 ， 主 葵 (Main 
Stem) 4S”; SEE (Tiller) WS, i 为 着 生 在 
EZE- RDE, 7 ATE BUT BE LAY 
TBAB. FHF MERE TEE RAB, 
本 研究 暂 不 考虑 。 首 先 在 主 葵 出 生 处 环绕 式 获取 
3 个 坐标 点 用 来 记录 茎 的 粗 度 信息 ， 然 后 从 茎 的 
出 生 处 开始 由 下 至 上 获取 若干 点 (图 3 (a) )， 
茎 上 最 后 一 片 叶 与 茎 的 连接 处 作为 结束 点 。 小 麦 
拔节 后 ， 要 求 每 个 点 与 小 麦 的 节点 位 置 重合 。 
穗 后 ， 需 以 穗 的 出 生 处 作为 结束 点 。 主 葵 数 据点 
集 由 {5 站 表示 ，N* 为 所 获取 主 茎 三 维 数字 化 点 
的 数量 ， 因 此 0<k<NY,ke2Z，2Z 为 整数 ; 将 
主 茎 出 生 点 SY 视 为 整个 小 麦 植株 的 出 生 点 。 分 
BEREA EMH SERR, NEII EAA 
据点 的 数量 。 

(2) 叶 。 由 工 (Leaf) KR, EZIK N 
L”; SEEI A L FEMA h E Abw tR 
方向 ， 逐 行 获取 叶片 左边 缘 、 叶 脉 、 右 边缘 三 个 
点 (图 3 (b) ) 数据 ， 以 叶 尖 处 一 个 点 作为 结 
束 。 要 求 间隔 距离 不 宜 过 长 ， 需 包含 能 明显 表示 
叶片 形态 变化 处 以 及 叶片 最 高 点 信息 ， 叶 片 发 生 
扭曲 时 仍 须 以 叶片 平展 开 之 后 的 初始 横向 方向 进 
行 数据 获取 。 主 茎 叶片 数据 点 集 由 {L”") 表 示 ， 
n 为 叶片 生长 发 育 顺 序 ，0, 为 主 芝 包含 叶片 数 
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Æ, KIEO <n< 06,n es Z; NA 为 所 获取 主 茎 第 
n 片 叶 三 维 数 字 化 点 的 数量 ,， 因 此 0<k< 
N", keZ; 分 蓝 时 片 数据 点 集 由 亿 人 人 表示 ， 
O, 为 对 应 分 昔 包 含 的 叶片 数量 ;，N2 ,为 对 应 分 
BESS n Fru lo KE A AE o 

(3) #8. HE (Ear) RI, FAWN E; 


(a) 茎 获取 的 数据 点 


分 右 的 穗 为 E“。 数 据 获取 从 穗 的 出 生 处 开始 由 
下 至 上 获取 若干 点 ， 以 穗 尖 处 作为 结束 点 。 主 蕉 
fA) OSCAR AUS rh (EE) SER, Nk PSII 
的 三 维 数字 化 点 的 数量 ， 因 此 
O<k<N™ ke Z; SBE BRAY BE Me h LEP) 


RI, NE ARI DLT BEA BE OE CEE o 


(b) 叶 片 获 取 的 数据 点 


图 3 小 麦 茎 叶片 数据 获取 示意 图 


Fig. 3 Data acquisition of wheat stems and leaves 


2.2 小 麦 表 型 参数 提取 方法 


2.2.1 数据 标准 化 

数据 可 视 化 和 提取 小 麦 植株 表 型 参数 前 ， 需 
对 三 维 数字 化 数据 进行 标准 化 处 理 ， 从 而 使 小 麦 
植株 均 位 于 原点 ， 且 主 茎 生长 方向 竖 直 问 上 重合 
于 Z 轴 ,便于 植株 的 可 视 化 显示 以 及 高 度 参数 的 
提取 。 

由 于 主 葵 数据 点 集 前 3 个 数据 点 为 粗 度 信 
息 ， 第 4 个 点 SY 为 主 茎 出 生 点 。 首 先 将 SY 平移 
至 坐标 〈0,0,0) ， 并 对 所 有 剩余 数据 点 做 同等 平 
移 操 作 。 然 后 根据 主 茎 生长 方向 对 植株 进行 旋 
转 ， 生 长 方向 由 主 共 上 第 一 个 和 最 后 一 个 记录 长 
度 信 息 的 坐标 点 形成 的 向 量 表 示 ， @ = SM - 
S¥z, v, = (0,0,1) 为 Z 轴 正方 向 上 一 向 量 ， 其 中 
Six, SM y, SM 2 Sb Bll Ze aN = BE ZS [ALS x. yz 坐标 
值 。 记 0 为 8 与 六 的 夹 角 ， 以 向 量 axis = gx v ME 
为 旋转 轴 ， 所 有 数据 点 绕 axis 旋 转角 度 0 完成 数 
据 标准 化 。 完 成 标准 化 之 后 ， 针 对 不 同 数 据点 集 
使 用 参数 提取 方法 进行 表 型 参数 的 计算 与 提取 。 


222 高度、 长 度 参 数 提取 

(1) 株 高 、 茎 长 、 节 间 长 的 计算 。 根 据 数据 
获取 规则 ， 拨 节 后 茎 的 数据 点 的 获取 与 节点 一 
BL. d (p,,p)) 表 示 三 维 空 间 中 两 点 的 欧式 距离 ， 
以 主 茎 为 例 ， 则 从 主 茎 出 生 处 开始 ， 第 i 个 节 间 
长 的 计算 方法 如 公式 (1) 所 示 。node 为 节 间 长 
BE, cmo 

node, = d(S".;,S"4), 1<i<N™-3 (1) 

葵 长 为 各 节 间 长 的 总 和 ， 计 算 方法 如 公式 (2) 
AN. sR, cm 

s= 2 node, (2) 

数据 标准 化 后 ， 三 维 空 间 内 的 Z 坐 标 就 表示 
高 度 信息 。 小 麦 植株 株 高 为 主 芭 和 分 慕 的 高 度 最 
大 值 ， 抽 穗 前 ， 高 度 最 大 值 由 最 上 部 叶 叶 尖 点 的 
Z 坐 标 决定 。N2"” 为 主 茎 最 上 部 叶 数 据点 的 数量 ， 
令 m = No", WW Lo ON Ee EBM A a 
坐标 点 ; N2。 Aaah Be T, AS ERRI AY ea 
坐标 点 ， 令 ;= Nilo. Wl Live’: oy Be ie EoD 
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叶 的 最 后 一 个 坐标 点 ， 株 高 的 计算 方法 如 公 
式 (3) Prax. HHR, em. 

H = max (Ly ®".z, max(L,"*"2)) (3) 

WR, tes Ip OK AAS EUT AS ee, ee 
PEMET AB TAY Z AB A EERIE o EM 
AEE EWS, Eye AEE T, E 
穗 的 最 后 一 个 坐标 点 ， 则 株 高 的 计算 方法 如 公 
式 (4) 所 示 。 


_ MO Typ Qr, M 
H = mar (L, "z max (L,” "zj Eyez, 


max (E wi .z)) (4) 

(2) 叶 长 、 叶 宽 、 叶 片 出 生 位 置 的 计算 。 数 
据 获 取 时 ， 叶 片 每 行 数据 点 的 间隔 距离 保证 了 曲 
率 没有 较 大 变化 ， 因 此 叶 长 采用 计算 叶脉 上 相 邻 
数据 点 间距 离 和 的 方式 。 但 折线 段 的 累积 和 小 于 
叶 长 实际 值 ， 根 据 不 同 品种 实际 值 和 提取 值 确定 
经 验 比 例 系 数 < 调 整 计 算 公式 。 由 于 集合 内 叶脉 
数据 点 的 间隔 为 3， 最 后 叶 尖 处 与 前 一 点 的 间隔 
为 2， 以 主 茎 第 一 片 叶 为 例 ， 数 据点 集合 大 小 为 
Ni ， 则 叶 长 的 计算 方法 如 公式 (5) 所 示 。 /为 


Ite, cm, 
Niw3 - 1 
at See 
i=l 
d (Lius Lyin) IC (5) 


叶 宽 的 计算 为 左 、 右 边缘 点 到 叶脉 点 的 距离 
之 和 ， 以 和 叶 长 同样 方式 确定 经 验 比例 系数 r 
以 主 茎 第 一 片 叶 为 例 ， 从 叶片 出 生 处 开始 ， 第 i 
排 位 置 处 叶 宽 的 计算 方法 如 公式 (6) 所 示 。w 
AM Si, cmo 

w; = (qd (Lala Lyi) + d (Lyn L3 ))/T, 

1 < i < N#”/3 (6) 

其 最 大 叶 宽 w 为 max (w,)。 

叶片 的 出 生 位 置 即 该 叶片 与 茎 的 连接 点 的 位 
置 ， 上 有 具体 为 第 1 行 数据 点 的 中 间 叶 脉 点 
(ig LA y, LY1z)， 且 出 生 高 度 为 L* zo 

(3) 穗 长 、 穗 出 生 位 置 的 计算 。 穗 长 也 采用 
穗 上 相 邻 数据 点 间距 离 和 的 方式 计算 ， 以 主 蔡 为 
例 ， 穗 长 的 计算 方法 如 公式 7) 所 示 ，e 为 穗 
K, cmo 


TE 
e= > d(E*,EM,) (7) 
i=l 


穗 的 出 生 位 置 为 穗 的 数据 点 中 的 第 一 个 坐标 点 ， 
有 具体 为 (EY.x, Ely, BY.z)， 且 出 生 高 度 为 EY.z。 
2.2.3” 粗 度 参数 提取 

(1) 茎 粗 的 计算 。 通 过 计算 茎 数据 点 集 前 3 
个 点 的 外 接 圆 半径 确定 茎 粗 。 根 据 三 点 共 面 、 三 
点 到 圆心 距离 相等 两 个 约束 条 件 求 得 3 个 平面 方 
程 如 (8) 所 示 。 

AP+B=0 (8) 

其 中 ，4 为 表示 方程 系数 的 三 阶 方 阵 ; BH 
常数 项 ; P 为 外 心 坐标 矩阵 。 以 主 茎 为 例 ， 通 过 
SXY、SN#、S4# 三 点 可 以 求 得 外 心 p,，(x,，y，2,)， 
则 蓉 粗 o 为 s*、SY#、SY# 三 点 到 圆心 距离 的 两 倍 
加 权 平 均值 ，cm， 计 算 方法 如 公式 (9) 所 示 。 


3 
o =£ Sa (5S, p.) (9) 
i=1 


(2) 围 度 的 计算 。 人 小麦 植株 的 紧凑 和 松散 程 
度 与 分 蓝 产 生 后 茎 秆 下 半 部 分 生长 情况 有 着 直接 
联系 ， 将 小 麦 下 部 1/3 处 各 葵 组 成 的 多 边 形 外 接 
圆 的 直径 定义 为 小 麦 的 围 度 。 株 高 为 五 的 小 麦 植 
株 通 过 拟 合 忌 3 tes BEND ESE mS Bop, FA 
于 坐标 点 不 在 确定 的 圆周 上 ， 所 以 利用 拉 格 朗 日 
乘 数 法 求解 圆心 的 最 优 解 ” 来 计算 小 麦 植株 的 
Eu (图 4)。 围 度 大 小 可 直接 说 明 小 麦 植株 的 
松散 程度 ， 定 量 描 述 小 麦 的 株 型 。 


HBr WA 


En 


图 4 植株 围 度 拟 合 示意 


Fig. 4 The fitting schematic diagram of wheat plant girth 


通过 构建 出 拉 格 明日 函数 已 (P,，4)， 分 别 
对 P,，4 求 一 阶 偏 导 数 ， 并 令 ”(P,,，4) 等 于 
Z, HERP = [x,y,z] " 即 为 最 后 圆心 坐标 的 
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最 优 解 ， 可 构建 方程 如 公式 (10) 所 示 。 
F(P,A)=|DP,- G| +A(UP, - 1) (10) 


其 中 ,DD 和 G 分 别 为 表示 线性 方程 组 f(x， 


yo z) TERAM OD x 3 阶 矩 阵 和 常数 


项 ; f ax yo z) 为 在 假设 所 有 坐标 点 都 在 圆 
周 上 此 前 提 下 ， 利 用 其 中 任意 两 点 的 连 线 垂 直 于 
连 线 中 点 与 圆心 连 线 此 性 质 所 构建 的 
MOD x 3 维 线性 方程 组 ; U= (abe) *, 
其 中 a，b，c 为 所 有 离散 点 (包括 圆心 ) 所 在 
HEH ax + by + cz- 1=0 的 方程 系数 。 

根据 各 茎 坐标 点 p, 到 圆心 p，(x,，y,，z,) 的 
两 倍加 权 平 均 距 离 确定 围 度 wx， 计算 方法 如 
公式 (11) 所 示 。 


2 n 
u =>) d(p»p.) (11) 
i=] 


2.2.4 角度 参数 提取 

(1) 茎 叶 夹 角 的 计算 。 根 据 茶叶 夹 角 定义 并 
参考 人 工 测 量 方法 ， 可 通过 叶片 出 生 处 生长 方向 
和 茎 局 部 生长 方向 这 两 个 空间 向 量 的 夹 角 计算 芸 
叶 夹 角 。 叶 片 出 生 处 生长 方向 1 可 由 前 2 个 叶脉 
数据 点 来 表示 ， 蔡 局 部 生长 方向 通过 叶片 出 生 高 
度 定义 此 叶片 在 第 几 节 间 ， 计 算 节 点 间 癌 量 s。 
则 茎 叶 夹 角 的 计算 如 公式 (12) tas, a HA 
Fa Cls 


Ur (12) 


a = arccos ( 
|| x 四 

(2) 叶 倾 角 的 计算 。 一 般 情况 下 ， 通 过 直接 
连接 叶 耳 与 叶 尖 ,测量 此 向 量 与 XOY 平 面 的 夹 
角 确 定 叶 倾角 。 但 小 麦 叶片 在 不 同时 期 存在 较 大 
的 差异 性 ， 弯 曲 程度 很 大 ， 所 以 需要 分 段 计算 。 
根据 数据 获取 规范 ， 叶 片 最 高 点 一 定 在 数据 点 集 
合 中 。 计 算 相 邻 叶脉 点 连 成 的 向 量 和 XOY 平 面 
之 间 的 夹 角 取 平 均值 ， 为 最 后 的 叶 倾角 结果 (图 5)。 
叶片 后 期 受 重力 原因 或 衰老 出 现下 垂 情况 ， 不 纳 
入 叶 倾角 计算 范围 内 。 以 主 葵 第 一 片 叶 为 例 ， 叶 
片 数据 点 集合 大 小 为 N， 所 以 (NM -=- 1)/3 Bui 


脉 向 量 中 只 取 Z 为 正 值 的 向 量 。 如 图 5 所 示 ， 叶 
片 最 高 点 之 前 的 有 效 分 段 向 量 (颜色 标记 段 ) 纳 
入 计算 范围 。 假 设 前 m 段 叶脉 向 量 符合 要 求 ， 取 
第 一 段 叶 脉 向 量 在 XOY 平 面 上 的 投影 向 量 @.， 
则 叶 倾 角 的 计算 如 公式 (13) 所 示 ，Z 为 叶 倾 
f, C) 


M, 1 M, 1 = 
(sa Zac-D+2) xe, 


) (13) 


1< 
= 一 > arccos 
B= = Drarecos 


M1 _ 7M,1 md 
| Lt, x |e| 


AS 小 麦 叶 片 叶 脉 分 段 向 量 示意 图 


Fig. 5 Segment vector and highest point of wheat leaf vein 


2.25 ”叶片 卷曲 .扭曲 程度 参数 提取 

在 小 麦 生长 过 程 中 ， 部 分 叶片 会 发 生 扭曲 、 
卷曲 等 形变 。 通 过 参数 对 叶片 扭曲 、 卷 曲 程度 进 
行 定量 化 ， 有 利于 分 析 不 同 品种 间 叶 片 的 形态 
差异 。 

(1) 叶片 卷曲 程度 的 计算 。 叶 面 卷曲 程度 刻 
画 以 叶脉 为 中 心 ， 左 右 两 侧 叶 面 形成 的 卷曲 。 通 
过 提取 叶 面 有 关 向 量 ， 定 义 小 麦 扭 曲 和 卷曲 程度 
参数 。 通 过 右手 定 则 计算 叶 面 数字 化 数据 中 叶脉 
AAIE they Al aig MESA y (图 6)， 从 而 确定 
此 叶 长 处 的 卷曲 程度 。 以 主 葵 第 一 片 叶 为 例 ， 叶 
片 数据 点 集合 大 小 为 Ni*， 则 能 确定 (NY - 1)/3 4 
叶脉 左右 面 元 夹 角 ， 卷 曲 程度 c 的 计算 如 公式 (14) 
所 示 。 

0 = BEES wo0 (14) 

c 取 值 范围 在 [0，1]j， 值 越 大 证 明 叶 片 卷曲 
越 明 显 。 

(2) 叶片 卷曲 程度 的 计算 。 叶 面 扭曲 程度 刻 
画 沿 叶脉 从 基部 到 叶 尖 方向 叶 面 方向 的 整体 变化 
趋势 。 未 发 生 扭曲 的 小 麦 叶 片 ， 对 应 的 左右 边缘 
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图 6 小 麦 叶 片 卷曲 程度 计算 示意 图 
Fig. 6 Schematic diagram of wheat leaf curling calculation 
点 确定 的 向 量 3 在 三 维 空间 中 应 该 具有 平行 的 性 
质 。 通 过 计算 相 邻 向 量 的 夹 角 累 积 和 (图 7)， 可 
以 确定 叶片 的 扭曲 程度 。 以 主 茎 第 一 片 叶 为 例 ， 
叶片 数据 点 集合 大 小 为 Nr*， 则 能 确定 (NY - 1)/3 
个 向 量 ， 向 量 计算 如 公式 (15) 所 示 。 


= 
= M, 1 M,1 
ve Ls go? 


1<ig(N” -1)/3 (15) 
扭曲 程度 gq 的 计算 方法 如 公式 (16) 所 示 。 
(NM - 1)/3 V, v 


q=( > arccos( ))/360 (16) 


$22 


图 7 小 麦 叶片 扭曲 程度 计算 示意 图 


Fig. 7 Schematic diagram of wheat leaf twisting calculation 


q 取 值 范围 在 [0，1]， 数 值 越 大 表明 叶片 扭 
曲 越 明显 。 图 8 给 出 了 对 小 麦 叶 片 数据 点 集 进行 
可 视 化 的 结果 ， 其 扭曲 度 分 别 为 0.12、0.56 
和 0.89。 


2.3 数据 可 视 化 


用 Open3D 库 对 小 麦 植株 数据 进行 处 理 并 可 
视 化 ， 数 据 处 理 流程 可 以 分 为 以 下 3 个 步骤 。 
(1) 三 角 网 格 化 。 利 用 数字 化 数据 包含 点 的 


/ 


| 
(a)0.12 (b)0.56 = (c)0.89 
图 8 三 种 扭曲 度 的 小 麦 叶片 


Fig. 8 Wheat leaves with three twisting degrees 


坐标 和 顺序 信息 ， 生 成 三 角 网 格 。 

(2) 网 格 细 分 。 对 三 角 网 格 化 后 的 叶片 ， 采 
H Loop 网 格 细 分 ”方法 对 各 叶片 网 格 加 细 ， 针 
对 小 麦 叶 片 沈 代 细 分 2 次 后 可 达到 最 佳 效 果 。 

(3) 网 格 着 色 。 根 据 小 麦 器 官 颜色 使 用 
paint_uniform_color 函数 对 网 格 进行 上 色 。 

图 9 分 别 给 出 了 小 麦 单个 叶片 和 整个 植株 原 
始 数字 化 数据 、 生 成 三 角 网 格 、 进 行 网 格 细 分 ， 
和 网 格 着 色 的 效果 图 。 


2.4 数据 分 析 方 法 


采用 以 下 几 个 标准 评价 参数 精度 。 

(1) 均 方 根 误差 (Root Mean Square Error, 
RMSE) 是 实测 值 与 预测 值 之 间 误 差 平 均值 的 平 
方 根 ， 计 算 方法 如 公式 (17) 所 示 。 


| S (actual, —predicted,)’ 
RMSE = |= 
A n 
(2) 决定 系数 (R-square, R) 是 实测 值 与 
预测 值 之 间 的 模型 拟 合 度 ， 取 值 范 围 在 10, 1], 


计算 方法 如 公式 (18) 所 示 。 


> (predicted, —Mean)’ 
二 | 


(17) 


R? (18) 


S (actual, -Mean } 
其 中 ，Mean 为 实测 值 均 值 ， 计 算 方 法 如 公 
式 (19) 所 示 。 


S actual, 
=I 


Mean = + 


(19) 
n 
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(a) 叶 片 原始 三 维 数字 化 数据 (b) 叶 片 三 角 网 格 化 


(e) 植 株 原 始 三 维 数字 化 数据 (人 植株 三 角 网 格 化 


(c) 叶 片 网 格 细 分 (dd) 叶片 网 格 着 色 


(h) 植 株 网 格 着 色 


(g) 植 株 网 格 细 分 


图 9 小 麦 叶 片 和 植株 可 视 化 数据 处 理 流 程 


Fig. 9 Visual data processing flow of wheat leaves and plants 


3 结果 与 分 析 


3.1 可 视 化 结果 


利用 Open3D 对 小 麦 三 维 数字 化 数据 进行 可 
视 化 ， 以 最 直观 的 方式 体现 小 麦 植 株 品 种 间 以 及 
不 同 生育 期 的 差异 ， 在 此 基础 上 进行 小 麦 模 型 构 
建 ， 可 解决 小 麦 絮 官 遮挡 严重 的 问题 。 

图 10 给 出 了 在 起 身 期 、 拔 节 期 、 抽 穗 期 获 
取 的 3 个 品种 数字 化 数据 的 可 视 化 结果 。 由 结果 
可 知 ，JM44 分 莹 最 多 ，XN979 次 之 ，FK13 最 
少 。 在 株 型 方面 ，XN979 在 拔节 前 植株 相对 更 
高 ; 但 拔节 后 FK13 节 间 迅速 伸 长 ， 直 至 形态 不 
再 发 生变 化 时 ，FK13 植株 均 为 最 高 ，XN979 相 
对 较 矮 ; 同时 FK13 围 度 值 较 大 ， 更 加 松散 ; 
XN979 和 JM44 围 度 值 较 小 ， 更 加 紧凑 。 在 叶片 
形态 方面 ，FK13 叶片 具 有 细 长 的 特点 ， 叶 长 最 
大 ， 叶 宽 最 小 ， 卷 曲 程度 较 大 ; XN979 和 JM44 
叶 宽 较 大 ， 卷 曲 程度 较 小 ;但 由 于 品种 特性 ， 
JM44 叶 片 拥有 更 大 的 叶片 扭曲 程度 ，XN979 叶 
片 最 为 平展 。 


3.2 ”参数 提取 结果 精度 分 析 


为 评价 所 提出 的 小 麦 植 株 表 型 参数 提取 方法 
的 精度 ， 选 取 了 起 身 期 、 拨 节 期 、 抽 穗 期 对 应 数 
据 进行 验证 。 其 他 方法 无 法 精准 提取 到 各 时 期 所 
需 表 型 参数 ， 因 此 利用 实测 数据 进行 对 比分 析 。 
表 1 给 出 了 抽穗 期 3 个 品种 小 麦 植株 长 度 、 粗 度 
和 角度 参数 的 提取 值 和 实测 值 均值 。 由 结果 可 
知 ， 提 取 值 与 实测 值 相差 不 大 ， 且 能 体现 出 品种 
差异 。 其 中 ，FK13 叶 长 值 最 大 ， 叶 宽 值 最 小 ; 
JM44 叶 片 拥有 相对 较 小 的 叶 长 值 ， 叶 宽 值 最 大 ; 
XN979 GH. AZK, AMM, FKI3AKA 
细 ，XN979 茎 偏 短 上 且 粗 ，JM44 居 中 。FK13 叶片 
的 荃 叶 夹 角 较 大 ， 但 叶 倾角 最 小 ; IM44 ZANT de 
角 最 小 ， 但 叶 倾 角 最 大 ; XN979 居 中 。 图 11 进 
一 步 给 出 了 3 个 品种 小 麦 植株 在 3 个 生育 期 提取 
的 叶 长 、 叶 宽 、 茎 长 、 茎 粗 、 茎 叶 夹 角 、 叶 倾角 
值 与 实测 值 对 比 结果 ， 以 及 验证 结果 均 方 根 误差 
RMSE 和 决定 系数 RR。 图 中 粉色 、 黑 色 和 蓝 色 分 
别 表示 起 身 期 、 拔 节 期 和 抽穗 期 三 个 生育 期 的 拟 
合 效果 。 
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© a 
a 丰 抗 13 号 (FK13) 西 农 979 号 (XN979) 济 麦 44 号 (JM44) 
CD 图 10 3 个 品种 小 麦 植株 (FK13、XN979、JM44) 在 不 同 生 育 期 (起 身 期 .拔节 期 抽穗 期 ) 数 据 可 视 化 结 
= Fig. 10 Data visualization results of 3 different cultivar wheat plants(FK13 .XN979 .JM44) in different growth periods (Rising 
N period, Jointing period, Heading period) 
> 表 1 3 个 品种 小 麦 植株 抽穗 期 时 长 .时 宽 、 茎 长 . 茎 粗 、 茎 在 小 麦 表 型 参数 提取 方面 ， 其 他 方法 都 无 法 
>< one eee 在 保持 植株 原 有 三 维 形态 结构 的 前 提 下 ， 进 行 小 
ie 1 _ ee and ia pi 麦 以 上 表 型 参数 的 完整 i 准确 获取 is, 特别 在 
width, stem length, stem diameter, angle between stem an _. ; rs 5 
— 三 叶 一 心 后 的 主要 关键 生育 时 期 ， 小 麦 形态 结构 
eal, and leat angie o iiferent cultivar whea ants 1n 
i ox leaf, and leaf angle of 3 different cult heat plant m 
(S) . | 复杂 ， 器 官 谈 挡 严 重 。 由 对 比 结果 可 知 ， 叶 长 、 
Heading period 
一 一 一 一 一 茎 长 、 茎 粗 、 茎 叶 夹 角 的 计算 误差 较 小 ， 三 个 品 
表 型 参数 均值 FK13 XN979 JM44 P >A mj 、 È 
提取 旗 叶 叶 长 /cm 15.07 14.56 14.26 种 的 平均 RR 分 别 为 0.93、0.98、0.93、0.85， 叶 
实测 旗 叶 叶 长 /cm 15.18 14.66 14.50 宽 和 叶 倾角 与 实测 值 的 相关 系数 尺 分 别 为 0.75、 
提取 旗 叶 叶 宽 /cm 1.11 1.82 2.06 0.73， 但 在 不 进行 植株 拆除 测量 的 情况 下 ， 误 差 
实测 旗 叶 叶 宽 /cm 1.13 1.79 1.96 不 大 。 
提取 茎 长 /cm 74.98 58.63 64.10 
er 4E no] ahr sz 
实测 茎 长 /cm 74.03 58.43 63.78 3.3 参数 提取 时 间 效 率 分 析 
提取 茎 粗 / 3.09 3.93 3.89 ses ` z 、 
ji 为 评价 本 方法 的 重建 效率 ， 计 算 了 单个 小 麦 
实测 茎 粗 /cm 3.11 3.90 3.92 植株 数据 处 理 和 参数 提取 的 CPU 所 用 时 间 ， 以 
ZB » TIA, 
提取 茎 叶 夹 角 /(°) 52.23 35.35 23.49 APRE | 
KMA - 夹 角 /(?) 53.86 35.21 23.40 抽穗 期 小 麦 为 例 ? 其 平均 分 莫 数 和 平均 叶片 数 分 
提取 叶 倾 角 /(°) 42.86 25.57 51.34 别 为 7 和 27， 在 此 情况 下 数据 读 取 和 人 处理 所 耗 时 间 
实测 叶 倾角 /人 ) 39.70 27.72 55.82 约 为 121 ms， 表 型 参数 提取 所 耗 时 间 约 为 12 ms, 


总 时 间 为 133 ms。 以 上 方法 是 在 windows10 操 作 
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图 11 3 个 品种 植株 在 不 同 生育 期 (起 身 期 拔节 期 抽穗 期 ) 提 取 的 叶 长 、 叶 宽 、 茎 长 、 
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Fig. 11 Comparison results between extracted and measured leaf length, leaf width, stem length, stem diameter, angle be- 


tween stem and leaf, and leaf angle of 3 different cultivar wheat plants in different growth periods( Rising, Jointing and 


系统 下 VS2019 中 实现 ， 使 用 的 计算 机 配置 为 : 4 
i7-10700 CPU @ 


Intel (R) (TM) 
2.90GHz, 16G RAM, 


Core 


Heading period ) 


讨 it 


ASTI FEET Y/N BE PUES ate EU PA BE 
征 ， 基 于 三 维 数字 化 技术 构建 了 一 套 快 速 、 精 
确 、 自 动 提取 小 麦 表 型 参数 的 方法 ， 具 有 一 致 性 
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强 、 精 度 高 、 操 作 便 捷 的 特点 。 此 方法 在 保持 植 
株 原 始 三 维 形 态 结 构 的 前 提 下 ， 对 小 麦 三 维 数字 
化 数据 进行 可 视 化 ， 并 一 次 性 计算 提取 小 麦 包含 
长 度 、 粗 度 和 角度 的 11 个 常规 可 测 表 型 参数 。 
同时 ， 定 量 分 析 小 麦 植株 的 松散 、 紧 凑 程 度 和 叶 
片 卷 曲 、 扭 曲 程度 ,分别 从 植株 和 器 官 等 不 同 尺 
度 更 好 地 描述 了 小 麦 的 形态 特征 。 此 方法 有 助 于 
分 析 小 麦 品种 的 形态 结构 差异 ， 对 小 麦 株 型 理想 
株 型 培育 具有 一 定 意 义 。 

利用 三 维 数字 化 仪 获取 单 株 小 麦 数 据 ， 获 取 
工作 量 与 对 所 有 参数 进行 人 工 测 量 相 近 。 但 利用 
三 维 数字 化 数据 可 对 如 植株 围 度 、 叶 片 卷曲 和 扭 
曲 程度 等 参数 定量 化 ， 人 工 则 难以 定量 测量 。 同 
时 ， 人 工 测 量 主观 性 强 ， 不 同 操作 、 个 人 读数 习 
惯 等 都 容易 影响 测量 结果 ， 而 三 维 数字 化 数据 更 
为 客观 ， 不 因 上 述 因素 产生 数据 误差 ， 具 有 更 好 
的 一 致 性 。 

基于 图 像 和 基于 三 维 点 云 的 植物 株 型 参数 提 
取 是 目前 应 用 较 多 的 两 类 方法 ， 可 实现 玉米 | 
棉花 、 甜 菜 ”等 植物 表 型 参数 的 提取 。 但 由 
于 小 麦 属于 多 分 莹 作物， 具有 分 蒙 和 叶片 多 、 叶 
片 细软 的 特点 ， 同 时 叶片 相互 遮挡 严重 ， 基 于 图 
像 和 点 云 的 方法 都 无 法 适用 于 小 麦 整 个 生育 期 ， 
且 提 取 的 表 型 参数 有 限 。 目 前 未 见 其 他 方法 在 不 
破坏 植株 原始 三 维 形态 结构 的 情况 下 ， 自 动 批量 
提取 本 系统 中 包含 的 所 有 小 麦 表 型 参数 。 在 已 有 
的 小 麦 表 型 参数 提取 研究 中 ， 利 用 多 视角 图 像 系 
统 采集 的 小 麦 幼 苗 图 片 ， 结 合 Mask R-CNN 方 
法 "提取 了 小 麦 叶 长 和 株 高 ， 相 关系 数 R 分 别 
为 0.87 和 0.98。 本 方法 提取 的 小 麦 叶 长 和 茎 长 RR 
分 别 为 0.93 和 0.98， 误 差 更 小 。 基 于 LiDAR 的 方 
法 "对 2016 一 2018 年 连续 两 年 小 麦 的 分 更 数 进 
行 了 估算 ， 相 关系 数 户 分 别 为 0.61 和 0.56， 同 时 
利用 高 度 指 标 和 体积 指标 构建 小 麦 生物 量 模型 ， 
但 缺乏 长 度 和 角度 等 更 精准 的 表 型 参数 的 获取 。 
本 方法 可 以 直接 根据 数据 集 数 量 得 知 分 蒙 数 ， 并 
提取 长 度 、 粗 度 、 角 度 相关 的 表 型 参数 。 提 取 的 
叶 长 、 茎 长 、 茎 粗 、 茎 叶 夹 角 的 尺 分 别 为 0.93、 


0.98, 0.93, 0.85, MH SE AN UF thi fA AY R Ad ll 
0.75、0.73。 其 中 ， 叶 宽 和 叶 倾角 的 误差 相 较 其 
他 表 型 参数 略 大 。 主 要 原因 是 小 麦 叶 片 较 窗 ， 同 
时 叶片 会 出 现 卷曲 和 扭曲 现象 ， 所 获取 的 叶 边 缘 
点 可 能 不 在 叶脉 的 垂 线 上 ， 人 工 测 量 时 也 不 易 确 
定 最 大 叶 宽 处 。 叶 倾角 误差 较 大 是 由 于 人 工 测量 
对 测量 位 置 选取 不 一 致 导致 测量 精度 不 高 ; 而 利 
用 三 维 数字 化 数据 计算 的 叶 倾角 是 多 处 提取 结果 
的 均值 ， 各 种 弯曲 程度 叶片 都 能 实现 精度 较 高 的 
计算 结 

本 方法 制定 了 包含 植株 语义 信息 的 数据 获取 
规范 ， 利 用 三 维 数字 化 仪 进行 数据 获取 精确 定位 
了 各 个 茎 、 叶 片 、 穗 等 器 官 的 空间 位 置 ， 解 决 了 
小 麦 器 官 间 遮挡 问题 ， 可 以 实现 小 麦 表 型 参数 的 
精准 提取 。 该 数据 获取 规范 可 直接 应 用 到 水 稳 等 
形态 结构 相似 的 禾 本 科 分 蒙 作物 上 ， 其 他 作物 也 
可 以 在 此 基础 上 根据 自身 形态 结构 特征 进行 修改 
和 应 用 。 

同时 ， 本 人 研究 实现 了 小 麦 植株 的 可 视 化 仿 
真 ， 能 完整 还 原 小 麦 三 维 形态 结构 并 进行 表 型 参 
数 的 精确 提取 。 虽 然 已 有 研究 实现 了 小 麦 三 维 模 
型 的 快速 重建 ， 但 由 于 利用 三 维 扫 描 或 多 视角 三 
维 重建 得 到 的 植株 点 云 存 在 如 叶 缘 和 植株 内 部 信 
息 缺 失 问题 ” ， 不 适用 于 复杂 植株 的 三 维 重 
建 ， 且 难以 高 精度 地 提取 器 官 尺 度 表 型 参 
数 光 "。 本 方法 获取 数据 精确 ， 但 依赖 于 仪器 设 
备 ， 数 据 获取 的 时 间 成 本 也 相对 较 高 ， 对 于 实现 
通 量 的 小 麦 表 型 参数 获取 还 需 进 一 步 研 究 。 
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Abstract: Aiming at the difficulty of accurately extract the phenotypic traits of plants and organs from images or point clouds 
caused by the multiple tillers and serious cross-occlusion among organs of wheat plants, to meet the needs of accurate phenotyp- 
ic analysis of wheat plants, three-dimensional (3D) digitization was used to extract phenotypic parameters of wheat plants. First- 
ly, digital representation method of wheat organs was given and a 3D digital data acquisition standard suitable for the whole 
growth period of wheat was formulated. According to this standard, data acquisition was carried out using a 3D digitizer. Based 
on the definition of phenotypic parameters and semantic coordinates information contained in the 3D digitizing data, eleven con- 
ventional measurable phenotypic parameters in three categories were quantitative extracted, including lengths, thicknesses, and 
angles of wheat plants and organs. Furthermore, two types of new parameters for shoot architecture and 3D leaf shape were de- 
fined. Plant girth was defined to quantitatively describe the looseness or compactness by fitting 3D discrete coordinates based 
on the least square method. For leaf shape, wheat leaf curling and twisting were defined and quantified according to the direc- 
tion change of leaf surface normal vector. Three wheat cultivars including FK13, XN979, and JM44 at three stages (rising stage, 
jointing stage, and heading stage) were used for method validation. The Open3D library was used to process and visualize 
wheat plant data. Visualization results showed that the acquired 3D digitization data of maize plants were realistic, and the data 
acquisition approach was capable to present morphological differences among different cultivars and growth stages. Validation 
results showed that the errors of stem length, leaf length, stem thickness, stem and leaf angle were relatively small. The R’ were 
0.93, 0.98, 0.93, and 0.85, respectively. The error of the leaf width and leaf inclination angle were also satisfactory, the R’ were 
0.75 and 0.73. Because wheat leaves are narrow and easy to curl, and some of the leaves have a large degree of bending, the er- 
ror of leaf width and leaf angle were relatively larger than other parameters. The data acquisition procedure was rather time-con- 
suming, while the data processing was quite efficient. It took around 133 ms to extract all mentioned parameters for a wheat 
plant containing 7 tillers and total 27 leaves. The proposed method could achieve convenient and accurate extraction of wheat 
phenotypes at individual plant and organ levels, and provide technical support for wheat shoot architecture related research. 

Key words: wheat; three-dimensional digitization; visualization; phenotypic traits extraction 
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